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ABSTRACT 

Two measurements of the cosmic-ray positron fraction as a function of energy have 
been made using the High Energy Antimatter Telescope (HEAT) balloon-borne instru- 
ment. The first flight took place from Ft. Sumner, New Mexico in 1994, and yielded 
results above the geomagnetic cutoff energy of 4.5 GeV. The second flight from Lynn 
Lake, Manitoba in 1995 permitted measurements over a larger energy interval, from 1 
GeV to 50 GeV. In this letter we present results on the positron fraction based on data 
from the Lynn Lake flight, and compare these with the previously published results from 
the Ft. Sumner flight. The results confirm that the positron fraction does not increase 
with energy above ~ 10 GeV, although a small excess above purely secondary produc- 
tion cannot be ruled out. At low energies the positron fraction is slightly larger than 
that reported from measurements made in the 1960 's. This effect could possibly be a 
consequence of charge dependence in the level of solar modulation. 

Subject headings: cosmic rays, elementary particles, instrumentatiomdetectors, ISM: 
abundances 
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1. Introduction 



suit. 



Cosmic-ray electrons and positrons interact with 
the interstellar medium exclusively through electro- 
magnetic processes, such as synchrotron radiation 
and inverse Compton scattering, which do not signif- 
icantly affect the nucleonic cosmic-ray components. 
For this reason, electrons are a unique probe of 
cosmic-ray confinement and source distribution in the 
galaxy. The observed flux is dominated by nega- 
tive electrons from primary acceleration sites. How- 
ever, about 10% of the total flux are secondary par- 
ticles, resulting from hadronic interactions between 
the nuclear cosmic rays and nuclei in the interstel- 
lar medium. These interactions produce electrons 
and positrons in roughly equal numbers. A num- 



ber of observations (Agrinier et al. 1969 




Buffing- 


;on et al. 1975 




Miillcr & Tang 1987 




Golden et al. 


1987, 


Golden et al. 1994 


) indicate that the positron 



fraction, e + /(e~ + e + ), increases with energy at en- 
ergies above 10 GeV. Such an increase would require 
either the appearance of a new source of positrons, or 
a depletion of primary electrons. Confirming cither 
of these possibilities would have a profound impact 
on our understanding of cosmic-ray sources. This 
motivated the construction of the HEAT (High En- 
ergy Antimatter Telescope) instrument, which was 
designed to determine the positron fraction over a 
wide energy range, and with improved statistical and 
systematic accuracy. A first balloon flight (HEAT- 
94) was launched in 1994 from Ft. Sumner, New 
Mexico, at a geomagnetic cutoff rigidity of roughly 



4.5 GV. The results from this flight (Bar wick et al 



1995 ) indicate that the positron fraction does not in- 
crease at high energies, and that positrons may well 
be of entirely secondary origin. A subsequent pub- 
lication (Golden et al. 1996) by another group re- 
ports a positron fraction which is statistically con- 
sistent with both the HEAT result and the previous 
measurements. 

In this letter we report the positron fraction mea- 
sured during the second HEAT flight (HEAT-95). 
This flight was launched in 1995 from a location with 
a low cutoff rigidity (Lynn Lake, Manitoba), in order 
to permit measurements over a larger energy range of 
1 GeV up to about 50 GeV. We present the results 
from this flight and compare this measurement with 
the previously reported data from HEAT-94. The two 
data sets have been analyzed in the same fashion, to 
avoid potential systematic differences in the final re- 



2. Detector Description 

The HEAT detector is described in detail elsewhere 
( Barwick et al. 1997 ). It has a geometrical accep- 



tance of 473 cm 2 -sr, and consists of a magnetic spec- 
trometer combined with a transition radiation detec- 
tor (TRD), an electromagnetic calorimeter (EC), and 
time-of-flight (ToF) scintillators. 

The magnetic spectrometer provides particle track- 
ing through an array of drift tubes (DT) in the mag- 
netic field generated by a two-coil, warm-bore super- 
conducting magnet. The DT array consists of 19 
tracking layers in the bending plane, and 8 layers in 
the non-bending plane, each providing a single-point 
tracking resolution of w 70/im. The performance of 
the magnetic spectrometer as a whole can be charac- 
terized by the maximum detectable rigidity (MDR), 
which is the rigidity at which the momentum of the 
particle is equal to the error in the momentum mea- 
surement. The mean MDR achieved in both flights is 
170 GV for electrons. 

While the magnetic spectrometer determines the 
momentum and the sign of the particle charge, ad- 
ditional measurements are needed to find the mag- 
nitude of the particle charge, the direction of parti- 
cle traversal through the instrument, the energy of 
electrons and positrons, and to reject hadronic back- 
ground. The TRD provides electron identification 
and hadronic background rejection through the de- 
tection of transition x-rays, which can only be pro- 
duced by particles possessing a large Lorentz factor 
(7 > 10 3 ). The TRD consists of 6 layers of propor- 
tional chambers and associated radiators. Using a 
maximum-likelihood technique we determine whether 
the signals induced on the TRD chamber cathodes 
indicate the presence of TR, or are consistent with 
ionization energy loss alone. In addition, the anode 
wires of each TRD chamber are read out in 25 ns time 
slices, in order to detect ionization clusters typical for 
x-ray signals. The time-slice data are analyzed with 
a neural-net technique and provide additional elec- 
tron discrimination. The combination of maximum 
likelihood and time-slice analysis achieves a hadron 
rejection factor of > 100 at an electron efficiency of 
90%. 

The EC consists of 10 layers of Lead and plas- 
tic scintillators. The 10 EC signals recorded for an 
event are used to measure the primary energy of the 
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particle, to determine the degree to which the lon- 
gitudinal shower development matches that expected 
for an electron, and the depth at which the shower 
starts. These quantities are obtained from a covari- 
ance analysis, using accelerator calibrations of the EC 



an d GEANT /FLUKA-based ( [Brun et al. 1994 l Fasso 
et frl. 1993] ) simulations. Over most of the energy 
range of interest the energy resolution of the EC is 
10%. 

The event trigger included the requirement that 
the energy deposited in the last 7 layers of the EC 
exceed that expected for a 0.5 GeV electron during 
the first flight, and a 1 GeV electron during the sec- 
ond flight. Although this requirement reduces the 
observed numbers of electrons and positrons near the 
threshold energy, it has no effect on the measured 
positron fraction, since the response of the EC to 
electrons and positrons of the same energy is iden- 
tical. The hadron rejection factor for the trigger and 
EC is 200 at an electron efficiency of 90%. An ad- 
ditional energy-dependent hadron rejection factor of 
1.5 to 10 is provided by a comparison between the en- 
ergy measured in the calorimeter with the momentum 
measured by the spectrometer. The hadron rejection 
factors and electron efficiencies of the EC and TRD 
are determined from the flight data by using the par- 
ticle identification obtained from one of the detectors 
to define clean samples of negative electrons and pro- 
tons, and applying the electron selection criteria of 
the other detector to these samples. 

The ToF system consists of a layer of four plastic 
scintillators at the top of the instrument, and the first 
3 scintillator layers of the EC. The top scintillators 
also provide a measurement of the magnitude of the 
particle charge to distinguish singly-charged particles 
from heavier nuclei. The time-of-flight measurement 
is used to eliminate upward-going (albedo) particles, 
which mimic antiparticles in the spectrometer. The 
rejection power of the ToF system is sufficient to re- 
duce the Helium and albedo background to negligible 
levels. 

3. Flight Summaries 

The first HEAT flight, HEAT-94, took place from 
Ft. Sumner, NM, on 1994 May 3. The total time 
at float altitude was 29.5 hours, with a mean atmo- 
spheric overburden of 5.7 g/cm 2 . During the flight, 
the geomagnetic cutoff rigidity at the detector varied 
between 4 GV and 4.5 GV. The second HEAT flight, 



HEAT-95, took place from Lynn Lake, Manitoba, on 
1995 August 23. During this flight, the total time at 
float altitude was 26 hours. The mean atmospheric 
overburden was 4.8 g/cm 2 , and the geomagnetic cut- 
off rigidity at the detector was well below 1 GV. The 
detector configuration was essentially the same for 
the two flights. In the HEAT-95 flight, one of the 
TRD chambers was inoperative due to a high- voltage 
system failure. The resulting loss of hadron rejec- 
tion power is small, and has been compensated for 
in the HEAT-95 analysis by slightly tighter electron 
selection criteria and correspondingly lower electron 
efficiency. All systems achieved comparable levels of 
performance in the two flights. 

4. Data analysis 

Both data sets are subjected to essentially identi- 
cal selection criteria to obtain a final sample of e^. 
A complete list of these selections is shown in Table 
1. The first set selects for singly- charged, downward- 
going particles which have a well-determined momen- 
tum, and a velocity measurement consistent with a 
(3 = 1 particle. The second set of criteria selects 
for e ± . The TRD maximum-likelihood and time-slice 
analyses, in addition to the shower-shape analysis ob- 
tained from the EC, result in a clean sample of e^. An 
energy and momentum selection appropriate to the 
geomagnetic rigidity cutoff of each flight is then made. 
Finally, we require that the energy, E , measured in 
the EC be consistent with the momentum, p , deter- 
mined with the magnetic spectrometer. We evaluate 
the measured distributions of the ratio E/p in order 
to determine the residual background in the data sets. 
The E/p ratio should be unity for e^, (subject to in- 
strumental resolution), but will normally have a value 
less than unity for hadrons. The E/p distributions ob- 
tained for the two flights are shown in Figure 1. The 
E/p interval used in the final data sets is shown as the 
cross-hatched region in this figure. The events falling 
outside this interval are primarily interacting hadrons 
which have survived the EC and TRD electron selec- 
tions, along with electrons and positrons falling in the 
low-side tail of their E/p distribution. Applying the 
E/p selection symmetrically ensures that the positron 
fraction is not biased by this selection. The level of 
background in the positron sample is estimated by de- 
termining the shape of the E/p distribution for inter- 
acting hadrons as well as that for electrons, and by fit- 
ting the measured E/p distribution for positron candi- 
dates to a superposition of these hadron and electron 
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distributions for each E/p interval. This background 
represents ss 1% of the positron counts at low ener- 
gies, increasing to almost 10% at high energies. 

Table 2 shows the corrected positron and electron 
counts and resulting positron fractions obtained from 
this analysis for the two data sets, binned according 
to the energy of the particle. The particle energy has 
been corrected for radiative losses to the top of the at- 
mosphere. The corrected electron and positron counts 
shown in Table 2 are obtained by subtracting the 
hadronic background and the secondary positrons and 
electrons produced in the atmosphere from the raw 
counts. The atmopheric contribution is determined 
by a Monte Carlo simulation of hadronic interactions 



of Aosmic rays in the atmosphere. Over most of the 
energy range of interest the flux of positrons and elec- 
trons produced in the atmopherc is found to be s=s 3% 
of the total e flux. An empirical estimate of the at- 
mospheric contribution is also obtained from the flight 
data by comparing the positron fraction measured at 
depths less than 4 g/cm 2 to that measured at depths 
greater than 6 g/cm 2 as a function of energy. The 
atmospheric corrections to the positron fraction de- 
termined from the Monte Carlo simulation and from 
direct measurements agree within the statistical error 
of the direct measurement. The systematic error in 
the positron fraction resulting from uncertainties in 
the atmospheric background correction is estimated 
to be 1% for energies well above the geomagnetic cut- 
off. 



5. Results 

The positron fraction as a function of energy is 
shown in Figure 2 for the two data sets, along with 
results from previous measurements by other groups, 
the predictions for purely secondary positron produc- 
tion, ( Prothcroe 82[ ), and the predicted positron frac- 
tion based on recent work ( Clem et al. 1996| ) which in- 
vestigated the possible effect of charge sign-dependent 
modulation. The HEAT-94 data points shown in Fig- 



ure 2 are essentially those previously published (Bar- 
wick et al. 1995). Minor changes, resulting from a re- 



finement of the atmospheric secondary correction, in- 
clude the elimination of events with energies between 
4.5 and 5 GeV, for which the secondary correction is 
more uncertain due to the proximity of the geomag- 
netic cutoff. The HEAT-95 measurement reinforces 
the conclusion that the positron fraction does not in- 
crease in the 10-50 GeV energy range: the results of 



the two measurements are consistent with each other, 
and with a general decrease of the positron fraction 
with energy. Both data sets do indicate an overabun- 
dance of positrons compared with the prediction of 
Prothcroe at all energies, but this disagreement may 
not be taken too seriously, as the model itself has 
inherent uncertainties. For example, the predicted 
positron fraction scales with the ratio of the absolute 
proton to electron flux, and the uncertainties in these 
quantities are directly reflected in the positron frac- 
tion. The combined data set suggests the presence of 
a feature in the positron fraction in the energy range 
from 7 to 20 GeV. Positro n production mechanisms 
have been suggested (e.g 



Kamionkowski & Turner 



1991 ) which would lead to an excess of positrons in 



this energy region, but the uncertainties in our data 
do not permit a definite conclusion, and further mea- 
surements are required to confirm that this feature 
exists. 

At energies below a few GeV, our measured posi- 
tron fraction is significantly higher than that reported 
in 1969 by Fanselow et al., while it is in excellent 
agreement with a recent measurement by an indepen- 
dent group ( Barbicllini et al. 1996| ). As the earlier 
measurement was performed at a different period in 
the solar activity cycle, it may be tempting to inter- 
pret the results as being affected by the charge depen- 
dence in the solar modulation. This effect, and its 
impact on the positron fraction measured at Earth, 
has been investigated by several groups (Moraal et 
al. 1991| , pern ct al. 1996| ). Clem et al. have de- 



veloped a model based on the observed systematic 
difference in the correlation between the electron flux 
measured in space by the ICE instrument and ground- 
based neutron monitor measurements for the 1980 
and 1990 solar polarity epochs. This model predicts 
that the present solar epoch favors the transmission of 
positive charge species, resulting in an enhancement 
in the measured positron fraction over the galactic 
value. While our observations may lend some sup- 
port to the hypothesis of charge sign-dependent solar 
modulation, a definitive test of this hypothesis will 
have to wait until the onset of the next solar epoch, 
about the year 2000. 
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Table 1 
Selection Criteria 



Selection Description 



HEAT-94 Selection Range HEAT-95 Selection Range 



TRD, DT track match 
Charge=l 
Velocity=c 
DT track x 2 
# track points 
DTH rigidity error 

TRD e ± M.L. 
# TRD chambers hit 
TRD time slice 
EC shower shape 
EC shower start 
Energy, momentum selection 
\E/p\ 



|Int DT -IntTHD !< 25 cm a 
0.77e < Z < 1.5e 
0.5< P <2.0 

x 2 < 10.0 

TO. fit > 9 

MDR/|R| > 4 

log(M.L) > 2 
Ntad =6 
Neural net output > 0.5 

xlc < 1-8 
X stari < 0.8 r.l. 
E > 3 GeV, | p |> 2.5 GeV/c 
0.7 < | E/p | < 3.0 



0.77e < Z < 1.5e 
0.8 < (3 < 2.0 
x 2 < 10-0 

n/it > 8 
MDR/|R| > 4 

log(M.L) > 2 
Ntad =5 
Neural net output > 0.9 
xlc < 1-8 

Xsiart < 0.5 r.l. 

E > 1 GeV, I p |> 1 GeV/c 
0.75 < | £7p I < 3.0 



a Int dt and IntTHD are the intercepts obtained in a fit of particle track in the DT and TRD 
systems 



Table 2 

± RESULTS 



Energy (GeV) (GeV) n» 



1.0 


- 1.5 


1.36 






65.9 


475.3 




0.122T0.016 


0.122T0.016 


1.5 


- 2.0 


1.76 






236.3 


1780.2 




0.117T0.008 


0.117T0.008 


2.0 


- 3.0 


2.46 






342.6 


3300.6 




0.094±0.005 


0.094T0.005 


3.0 


- 4.5 


3.64 






205.5 


2631.6 




0.072T0.005 


0.072T0.005 


4.5 


- 6.0 


5.14 


48.6 


730.5 


62.6 


1218.6 


0.062T0.010 


0.049±0.006 


0.054T0.006 


6.0 


- 8.9 


7.20 


90.1 


1049.3 


48.8 


846.5 


0.079±0.008 


0.055T0.008 


0.068T0.006 


8.9 - 


- 14.8 


11.0 


38.6 


571.9 


20.0 


455.6 


0.063±0.010 


0.060±0.012 


0.062±0.008 


14.8 


- 26.5 


18.4 


13.7 


227.8 


6.9 


148.1 


o.o57±8:8il 


o.o44±8:8?l 


0.052±0.013 


26.5 


- 50.0 


32.3 


2.1 


41.2 


2.1 


29.4 


0.048±m? 


0.070±m 


0.057±8;8P 



1 / = e+/(e-+e+) 
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Fig. 1. — E/p distributions for e ± selection. The dashed curve is a Monte Carlo calculation of the expected 
distributions. The cross-hatched region indicates the E/p interval used in the final data sample. Note that the 
vertical scale is logarithmic. 
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Fig. 2. — The positron fraction as a function of energy obtained from the two data sets. The error bars applied to 
the HEAT data points represent statistical errors only. 
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